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Abstract 

Liouvillian systems were initially introduced in [3] and can be seen as a 
' natural extension of differential flat systems. Many physical non flat systems 

. seem to be Liouvillian (cf [12j [ITJ Q] [10] ) . We present in this paper an alter- 

^ ' native definition to this class of systems using the language of diffieties and 

O . infinite prolongation theory. 

1 Introduction 

OV 

Liouvillian systems were initially denned in the differential algebra setting. We 
give here a new formulation using the language of diffieties and infinite dimensional 
geometries. This mathematical framework is well suited to study Liouvillian sys- 
tems. Recall that one of the main property of flat systems is that the variables 
of the system (state, inputs) can be directly expressed, without any integration 
of differential equations, in terms of the flat output and a finite number of its 
time derivative. Liouvillian systems share a similar property. To be able to derive 
the trajectories of a Liouvillian system, we also need some elementary integrations 
called quadratures. This can be illustrate through the following academic example 

X2 = X 3 , (1) 
X 3 = U. 

It is quite easy to show that (P) is flat for i = 1,2 and not flat for i = 3 ([2])- 
However, for i = 3, the subsystem 



CO 

d 



X2 = X 3 , 
X 3 = U. 



(2) 



*The author would like to thank Prof. Michel Fliess for helpful discussions and his constant 
support during this work. 
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is flat with a flat output y = %2 and the trajectory of x\ can be obtained by mean 
of an elementary integration x\ = J y + y 2 . 

For the sake of convenience, we first recall, in sections [2] and El some facts 
concerning the theory of diffieties and the Lie-Backlund approach to equivalence 
and flatness (cf [H El El 13 E3] ) • In section HI we define Liouvillian systems using 
the language of diffieties. Finally, we illustrate the class of Liouvillian through the 
concrete case of rolling bodies (pp, [10], [2]). 

2 The language of diffieties 

Let / be a countable set of cardinality £, which may be finite or not, and the 
linear space of all real-valued functions x = (x % ) on /. The space MJ has the 
natural topology of the Euclidean space if / is finite and the Frechet topology 
otherwise. The elements x l , i £ I, are called coordinates. For an open set U C M. 1 
we denote by C°°(U) the space of all real- valued functions on U that depend 
on finitely many coordinates and are smooth as functions of a finite number of 
variables. A chart on a set M is a 3-tuple (E7,<p,R 7 ), where U is a subset of M, 
ip is a bijection of U onto an open subset <p(U). The notions of smooth charts 
and smooth atlases can be defined as in the finite dimensional case. The set M, 
equipped with an equivalence class of smooth atlases, is called a C°° M 7 -manifold. 
The number £ does not depend on a chart (U, ¥?,M 7 ) and is called the dimension 
of the smooth manifold M. 

A diffiety is a pair Jt = (M, CTM) where M is a C°° M 7 -manifold and CTM 
a finite dimensional involutive distribution on M. The distribution CTM is called 
Cartan distribution and its dimension the Cartan dimension of j% . Local smooth 
sections of CTM are called Cartan fields. We are only concerned here with the 
case of ordinary diffieties, i.e., the dimension of CTM is equal to 1. For the 
sake of convenience, we use without distinction the notations (M, CTM) and 
(M, 8m) to denote the ordinary diffiety where 8m is a basis vector field of 
CTM. Let Jt = (M, CTM) be a diffiety with dim CTM = 1. Let (C/,^,^) 
be a chart on M and 8m be a basis vector field of CTM on U, then the 4-tuple 
(U, if, M 7 , 8m) is called a chart on We denote by kert?M the kernel of the 
linear map 8 M ■ C°°{M) -> C°°(M), i.e., 

ker8 M = {$ € C°°(M) / 8 M <& = 0}. 

A real- valued C°° function i9 on M such that i? G ker 8m is called a local first 
integral on Ad. A local first integral is said to be trivial if it is a constant ([§]). 
Let : M —> N be a smooth mapping. As usual, we denote by 0* : TM — >• TiV 
the differential (or tangent) mapping of (f>, where TM (resp. TN) is the tangent 



bundle of M (resp. N), and by (ft* : T*N — > T*M the dual differential mapping 
of (ft, i.e., the dual mapping of (ft*, where T*M (resp. T*N) cotangent bundle of 
M (resp. N). 

A smooth mapping (ft : M — >■ N is called a Lie-Backlund morphism of a diffiety 
Jl = (M, CTM) into a diffiety ^ = (JV, CTN), written (ft : ^ JV , if it 
is compatible with the Cartan distributions CTM and CTN , i.e., (ft*(CTM) C 
CTJV. 

Example 2.1. Denote by R™ = R m x R m x • • • £/ie product of a countably infinite 
number of copies o/R m . Consider the ordinary diffiety & = (F,CTF), where 
F = lx R™, and let (U,ip,M. x RJ£,c?f) be a chart on & with local coordinates 
{t, w\ | i = 1, . . . , m ; ^ > 0} and frasis Cartan field 



am o 



The diffiety J^", as above defined, is usually called trivial diffiety and plays a central 
role in the Lie-Backlund approach of flatness. On some occasions, we will use the 
short notation 

to represent the basis Cartan field 8f- ■ 



3 Flat systems 

A diffiety is said to be (locally) of finite type if there exists a (local) Lie- 
Backlund submersion it : ^ — > & such that the fibers are finite dimensional. The 
integer m is called the (local) differential dimension of ^# (cf. [5]). 

Definition 3.1 ([5j [6]). A system is a (local) Lie-Backlund fiber bundle om = 
(^#, R, A), where 

(i) jM is a diffiety of finite type where a Cartan field &m has been chosen once 
for all; 

(ii) M is endowed with a canonical structure of a diffiety, with global coordinate t 
and Cartan field d/dt; 

(Hi) A : ^# — > R is a Lie-Backlund submersion such that \*(8m) = d/dt. 



□ 



The system of = R, pr), where pr is the natural projection mapping 
pr : {t,w^} t and & a trivial diffiety, is called a trivial system. 

The differential dimension of a system om = (^#,K, A), denoted dim cliff <tm, 
is the differential dimension of the associated dimety ■ 

Definition 3.2 (|5J). Two systems a M = (^#,M, A) and <jn = (^V,R,d) are said 
to be (differentially) equivalent, written o~m — 0~n, if and only if 

(i) <^>*(Sm) = 9n, where 4> ■ ^ ^ is a Lie-Backlund isomorphism; 

(ii) A = 4>*5. 

□ 

A system om = (^#, K, A) is said to be (locally) differentially flat, or simply 
flat if it is (locally) equivalent to a trivial system. If {t, yf | i = 1, . . . , m ; ^>0} 
are local coordinates on & then y = (yi, ■ ■ ■ ,y m ) is called a flat or linearizing 
output. 

4 Liouvillian Systems 

A diffiety S* = {S, CTS) is called a subdiffiety of a diffiety Jt = (M, CTM) if 

5 is a submanifold of M and CTS = TS n CT^M, i.e., the natural embedding 
i : J^ 7 — > ^ is a Lie-Backlund immersion^. The fiber bundle T$M denotes here 
the restriction of the vector bundle TM on S, i.e., 

T S M = |J T p M. 

P&S 

The tangent mapping : TS — > TM is injective and the image l*(TS) C T$M. 
If ^ is of finite type, then clearly 5? is of finite type as well. 

Definition 4.1. ^4 system cfm = (^,M, A) is said to be a differential extension 
of as = («5^,R, <5), denoted as/o~M or as C (Jm, onZy i/ 

(i) y is a subdiffiety of '; 

(ii) the restriction t*A = 5, where i* is the dual mapping of the natural embedding 



1 Since we consider only diffieties of Cartan dimension 1, CTS = CTsM here. 



□ 



Consider the differential extension om/°~S-, with dim differ m = m - Since 
is of finite type, there exists a Lie-Backlund submersion ir : d% — > & such that 
its fibers are finite dimensional, say re. Assume now that au is not flat and as is 
a flat with a flat output y = (yi, . . . , y m ). Define the canonical bundle morphism 
p : TM — > TM/TS that takes a vector £ € T p M, p € M, to its equivalence 
class ( + T p S and let r : TM/TS — > M be the fiber bundle whose fibers r _1 (p), 
p € M, are finite dimensional. If {t, 771, ... , r? s , | i = 1, . . . , m ; v > 0} are local 
coordinates on 5? then the Cartan distribution of 5? is spanned by 

ds o m o 

^ „i 9 x-x- (y+i) a 



* fe£S 4 a^' 



.7 

where iy are C°° functions on S. Using the short notation, ds can be written 
under the form 



at dri a^w 



with F 1 = (F^, . . . , F^). A local smooth section £ of TM/TS is given by 



^ j ae, ac 



c= V F 2 — = F 2 



where F 2 are C°° functions on M, F 2 = (F 2 , . . . , Fj), with 

••• >€d,m,-- -,Vs,Ui \i = l,...,m; v > 0} 

local coordinates on 

Definition 4.2. Let ctm 5e a differential extension of a flat system as and y a 
given flat output of as- Then, as is called a flat subsystem of oyi and the flat 
output y of as, a partial linearizing output ofau- If, in addition, that flat output y 
is such that d = dimr _1 (p), with p € M and r : TM/TS — > M the aforementioned 
fiber bundle, is minimal, then d is called the defect, as a maximal flat subsystem 
and y a maximal linearizing output of ' o~m- D 

Consider now the classical dynamics 

x = F{x, u), (x, u) e X x U C R n x R m , (3) 

where x = (x±, . . . , x n ), u = {u\, . . . , u m ) and F = (Fi, . . . , F n ) is a m-tuple of C°° 
functions on XxU. To © we can associate a diffiety Jt = (Rxl n xM m xl™, d M ) 
with local coordinates {t, x%, . . . , x n: uf^ \ i = 1, . . . , m; v > 0} and Cartan field 



A subsystem of (J3J) is given by a diffiety 5? = (S,ds), with local coordinates 

{i,77i, . . .,r] s ,uP \ i = l,...,m; v > 0} 
and a basis Cartan field 

at j=i ar? J i=l^>0 <9u> 

where r\ = (rji, . . . , n s ) £ X 1 C R s and are C°° functions on X 1 x U. A local 
section £ of TM/TS is given by 

d=n—s n q 

C= E if(^,u)^- = F 2 (*,«)-, 

where £ = . . . € X 2 C M d , F 2 = (Ff, . . . and F 2 are C°° functions 
on X 1 x X 2 x [/ = I x [/, The vector £ represents only the complement of 77 
(by renumbering the Xj's if needed) to form the vector x, i.e., x = (?7,£)- We 
can assume, in the sequel, that coordinates rj and £ are given by the projection 
mappings pv 1 and pr 2 

pri : {t,x 1 ,...,x n ,u^} -)• {t,r?i,...,r? s ,nf ) } 
pr 2 : {t,xx,...,x n ,uf>} ->• {*, 6, •••,&} 

Thus, if o"m is a differential extension of a flat system as, then dynamics (|3|) 
admits the following decomposition 



x 



Definition 4.3. A system om is called a Piccard-Vessiot extension of a system 
OS if • 

(i) as is flat; 

(ii) the Cartan field 8m is of the form 

d 

d M = A{r),u)£— +d s . 

(in) kerdM = kerdg. 




A differential extension omI°~s such that om is a Piccard-Vessiot extension of as 
is said to be a Piccard-Vessiot extension. □ 

In the case of a Piccard-Vessiot system, (|3|) takes the form 

(4) 




where A(rj,u) is a n x n matrix of C°°(X 1 x U) functions. 

Proposition 4.1. LetaM/&s be a Piccard-Vessiot extension. Then, au is locally 
controllable. □ 

Proof. Since as is flat, as — ap, hence ker (9^ = kerds = R, i.e., any local first 
integral of au is trivial, and it follows that au is locally controllable (cf [8]). □ 

Definition 4.4. Let au be a differential extension of a flat system as and y 
a flat output of as- Then, au is said to be a Liouvillian extension of as, or 
simply aul^s is a Liouvillian extension, if and only if there exists a nested chain 
of subsystems as = as C <7Si C ••• C as d = au, with as ] = {S?j,M.,8j) and 
S^j = (Sj , dsj ), such that, for j = 1, . . . , d, ker dsj = ker ds j _ 1 , where either 



(i) 8 Sj = ajd/dSj + d 3j _ lt a, € C°°(S3-i), or 

(a) d Sj = ojtjO OZj ■ d Sj ,. a: e c°°(Sj ,;. 

If as is maximal (resp. partial), i.e., d is the defect of au, then au is called 
Liouvillian system (resp. partial Liouvillian system) and y Liouvillian output 
(resp. partial Liouvillian output). □ 

Remark 4.1. According to the definition, a local section (j ofTSj/TSj-i is given 
either by 

00 Cj = Uj(r),£i,...,£ j -i,u)d/d€ j (hence £j = aj(r], fr, . . . , u) and ^ = 
(ii) (j = aj(r],Ci, ■ ■ ■ ,^j-i,u)Cjd/d^j (hence £j = 0^(77, . . . , ^_a, u)£j and ^ = 

Hence, Liouvillian extensions are extensions by integrals (i) or exponential of in- 
tegral (ii), usually called extensions by quadratures. ■ 



Actually, it is easy to see that Liouvillian extensions are a particular case 
of Piccard-Vessiot extensions. For (ii) of definition 14.44 it 1S clear that asj is a 
Piccard-Vessiot extension of o~s j „ 1 - The extension by integral can be obtained by 
considering the Piccard-Vessiot extension agj of (?Sj-i with Cartan field given by 

with etj = dctj/dt and atj G C°°(Sj). 

Remark 4.2. Notice that an arbitrary linearizing output y for as does not neces- 
sarily give rise to a Liouvillian system. Therefore, the Liouvillian character of a 
system strongly depends on the choice of y. ■ 

Let T n {C°° (X 1 x W 71 )) the set of n x n lower triangular matrices with compo- 
nents in C°°{X l x R m ) and U^C 00 ^ 1 x R m )) the subset of T n (C°°{X l x R m )) 
such that all the diagonal components are equals to 1. 

Theorem 4.1. Consider the Piccard-Vessiot system Q). IfA(n,u) G U<i(C °(X 1 x 
W 71 )) then is Liouvillian. □ 

Proof. Let A = (a i): ,) iji=M G U^C^iX 1 X R m )), then 

i 

6 = X a *.i^i' witn °m = 1 ' * = 1 ' • • • ' d - 

5=1 

In particular, 

6 = 6, 

and it follows that 6 is an exponential of integral. Next, 

1?\ 6 6 
a J 6 a 



in the other words, 

6 = 6 y «2,i- 

Finally, differentiating i = 1, . . . , d, gives 



6~\ . 6 . . 6-i 

7- ) = + «i,27- + • • • + <H,i-V 



6 y 6 ? 



i 



Now, making the appropriate induction assumption, we deduce that £j, i 
2, . . . , d, can be obtained by mean of the integral 

t t [ , & , i &-i 

Si = Si / + a i,2"^ + . . . + GUi-1— — , 

J SI 51 

which concludes the proof. □ 

Theorem 4.2. Consider the Piccard-Vessiot system If A(rj, u) £ Td(C°°(X 1 
R m )) then is Liouvillian. 

Proof. Let A = (a id ) hJ=l4 6 T d (C°°{X x x R m )), i.e., 

i 

ii = ^a itj ij, i = l,...,d. 
3=1 

and set 

Ci — Oi,i, i — lj • • • j d. 

So the q's are integrals of elements of the flat subsystem as- For i = 1, we get 

Cl = ai,i£l> 

and it follows that £i = e Cl . Next, 

M 6 6 6 . , \ 6 

7- ) = J - - 7TF = ° 2 . 1 + ' a2 > 2 ~ a ^T' 
SI y SI SI SI SI 

in the other words, 

6 = 6e C2 - Cl |a 2il e Cl - C2 . 
Finally, differentiating i = 2, . . . , d, gives 

, S*2 . . . / 

7- I = Oi,i + a i)2 — + • • • + ai,i_i— h (aj,, - ai.ij — • 

Si y Si Si Si 

Now, making the appropriate induction assumption, we deduce that £j, i 
2, . . . , d, can be obtained by mean of the relation 

ii = ti^~ Cl • J (o<,i + a i>2 | + . . . + Oi.i-i^i) e c i- c S 

which concludes the proof. □ 



5 The rolling bodies 



Let us illustrate the class of Liouvillian through the concrete case of rolling bodies 
(see [1] for a more larger treatment). The kinematic equations of motion of the 
contact point between two bodies rolling on top of each other are given in geodesic 
coordinates by 

vi = Ul, 

V2 = ui cos ip — U2 simp , (5) 
W2 = — j?(iti sm.ip + U2 cosip), 
^ = -ff^U2 — ^ff-(ui simp + U2 cos 

where B = B(vi,wi), C = C{v2,W2), B Vl = dB/dvi and C V2 = dC/dv 2 . In the 
case of the well known plate ball problem, B = 1 and C = cos{v 2 ) and ([5]) takes 
the form 

vi = Ul, 
Wl = u 2 , 

i)2 = ui cos i\) — U2 sin ijj, (6) 
^2 = -■^^[(.uisinif) + u 2 cosip), 
ip = tan V2 (ui sin ijj + U2 cos ip) . 

This system is Liouvillian and a Liouvillian output is given by (see pQ for the 
details) 

x = vi-v 2 cosip, 
y = wi + w 2 sinip. 

The idea to give a new formulation of Liouvillian systems within the mathematical 
framework of dimeties was actually motivated by the study of the rolling bodies 
system. As matter of fact, it is not clear whether this system is Liouvillian within 
the differential algebraic setting (cf [3]). 

First, notice that system [6] is not under a suitable form to describe a system 
in the differential algebraic setting (the associated differential field extension in 
not finitely generated). However, using transformations a = tan(^/2) and £ = 
tan (v 2/2), ([5]) writes (plate ball case) 

vi = Ul, 

Wl = U 2 , 

Z = W&) [C 1 " o 2 )*! ~ 2au A . (8) 



W2 = ~ { i_^i +a 2) [2CTU1 + (1 ~ <T 2 )U2 

& = jS^[2aui + (1 - a 2 )u 2 



and becomes explicit and rational. The associated differential ideal is thus prime 
and leads to finitely generated differential field extension. Nevertheless, the Liou- 
villian output ([7]) takes now the form 

x = v\ — 2arctan£, 

y = w 1 + j^jw 2 . 

If we denote by M.(v i, wi, £, W2, cr, u\, u 2 ) the differential field generated by R and 
the variables {v\,w\,^,W2,o~, U\,U2\, and w%, £, W2, cr, u\, 112) its algebraic 

closure, then x and y are not in u>x, £, u>2, cr, u%, U2), and it follows that Q 
is not a Liouvillian output for ([5J) in this context. 



0) 
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